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Abstract 

In this paper we discuss some topics related to the general theory of frames. In 
particular we focus our attention to the existence of different 'reconstruction for- 
mulas' for a given vector of a certain Hilbert space and to some refinement of the 
perturbative approach for the computation of the dual frame. 



1 Introduction 



Whenever we deal with a (separable) Hilbert space TC the first problem we face with is 
the way in which an arbitrary element / ETC can be conveniently expressed. As we know 
the usual choice is to expand / in terms of an orthonormal basis {e n } of TC: in this way 
the expansion is particularly simple, / = J2n < e n, f > e n , and the Parseval equality is 
satisfied, J2n I < e n, f > | 2 = ll/ll- It is also well known that in any Hilbert space we can 
build up different orthonormal bases, for instance using unitary transformations starting 
with one of such basis. 

Sometimes, however, the conditions of our (physical) problem force us to consider a set 
of vectors which is no longer orthonormal but is still a basis of TC. We have a typical 
example of this situation when the set forms what is known as Riesz basis, see 

and references therein. One possible characterization of a Riesz basis is the following: 

The set is a Riesz basis if and only if 

- all the $ n are independent, that is no $„ o lies within the closure of the finite linear 
span of the other $ n , and 

- 3 A > 0, £> < oo so that, for any f E TC, 

A||/|| 2 <EI<^,/>| 2 <5||/|| 2 . (i.i) 

n 

This last property implies that the vectors of the set {$ n } generates the whole Hilbert 
space, while the first condition says that these vectors are linearly independent: therefore 
they form a typical basis of TC. 

We know many examples of sets of vectors which are not a Riesz basis but still have a 
relevant role in the description of vectors of TC. For instance, if we consider any overcom- 
plete set of coherent states, f2|, this set satisfies a relation similar to the one in (11.1ft , but 
the vectors are not linearly independent. Equation ( |1 . 1| ) is also satisfied by some sets of 
wavelets, . Sets of vectors of this kind are known as frames. In other words, we can say 
that a frame is a set of generators of TC, but, since the vectors are not independent, the 
way in which a vector / ETC can be expanded by means of these vectors is not unique. 

Many properties of the frames have been discussed in the literature. In particular the 
interest of the researchers has been focused mainly on the reconstruction formula: given 
a frame {<p n } in which way the sequence {< <p n , f >} is linked to the function /? What 
does the Parseval equality become for frames? It is possible to operate in such a way to 
optimize the numerical procedure? 

All these questions have already an answer in the literature. In this paper we will 
show that many other relevant questions are still to be raised and can be solved: it is 
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possible to find different reconstruction formulas, maybe more convenient than the usual 
one? It is possible to build up a perturbation scheme which works in conditions in which 
the usual one, @, has difficulties? 
The paper is organized as follows: 

in Section 2 we introduce a bit of notation and, in order to keep the paper self con- 
tained, we quickly recall the main known results about frames. 

In Section 3 we show how, starting from a given frame, it is always possible to define 
a tight frame with frame bound equal to one, as well as many other frames which give 
different reconstruction formulas. We will discuss many examples of this procedure and 
we will also show how the usual perturbation scheme works in this case. 

In Section 4 we discuss a different perturbation approach which may be useful whenever 
the frame constants are not close to each other. 

In Section 5 we discuss some analogies between frames and orthonormal bases. 

2 Notation and Known Results 

Most of the results we are going to summarize in this Section can be found in || and |4| . 
We also like to cite || which is an useful review on wavelets with a long Section on frames 
and where it is contained an excellent bibliography. 

Let "H be an Hilbert space and J a given set of indexes. Let also X = {(f n (x)}, n G J, 
be a set of vectors of 7i. We say that X is an (A,B)-/rame of 7i if there exist two positive 
constants, called frame bounds, A > and B < oo, such that the inequalities 

^imi 2 <Ei<^/>i 2 < 5 imi 2 (2- 1 ) 

n 

hold for any f E7i. 

To any such set X can be associated a bounded operator F : 7i — > 1 2 (J) = {{c„}„ e j : 
J2n&j \ c n\ 2 } < oo defined by the formula 

VfEH {Ff) j =<cp j J>. (2.2) 

We will omit the dependence on X of F whenever this does not cause confusion. Due to 
equation fl2.1p we see that \\F\\ < \[B. The adjoint of the operator F, F*, which maps 
1 2 (J) into TC, is such that 

V{c}Gl 2 (J) F*c = Y,*Vi ( 2 -3) 
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Of course, also F* depends on the family X we are considering. Again, we will omit this 
dependence whenever it will not cause confusion. 

By means of these operators condition ( |2.1|) can be rewritten in the following equivalent 
way: X is an (A,B)-/rame of TC if there exist two positive constants, A > and B < oo, 
such that the inequalities 

Al < F*F < Bl (2.4) 

hold in the sense of the operators, |§. We have used 1 to identify the identity operator 
in B(H). 

Condition fl2.4|) implies that the operator F*F, which maps Tt into itself, can be 
inverted and that its inverse, (F*F)~ 1 is still bounded in TC. In other terms, we have that 
both F*F and (F*F)- 1 belong to B(TC). 

Following the literature, see || for instance, one define the dual frame of X, X, the set 
of vectors ifii defined by 

(p. = (F*F)- 1 i Pl Vt G J. (2.5) 

We have called X the dual frame since, as a matter of fact, it is really a frame, and, in 
particular, it is a (-^, ^)-frame, see Defining now a new operator between TC and 1 2 (J) 
as F = F(F*F)~ 1 , it is easy to prove that F is such that (Ff)i =< fy, f >, for all f ETC. 
Moreover, the following relations hold: F*F = F*F = 1. By means of these equalities it 
is proved that any vector of the Hilbert space can be expanded as linear combinations of 
the vectors of the set X or of the set X. We have the following reconstruction formulas: 

/ = < Vuf > & = < <Piif > Vi ( 2 -6) 
ieJ ieJ 

for all f £ TC. In it is also discussed that, since ifi = (pi for all i G J, then the dual 
frame of the set X is nothing but the set X itself. 

As it is clear from equation (|2.6| ), a crucial role in the reconstruction procedure is the 
knowledge of the set X. In order to obtain the explicit expression for (p^, we first have to 
know how the operator [F*F)~ l acts on the vectors of TC. This is, in general, a difficult 
problem to solve. Only in a single situation we can give an easy answer, namely when 
our frame is tight. This means that the frame bounds A and B coincide, A = B, so that 
equation (|2.4| ) reduces to 

F*F = AH, (2.7) 

which implies also that (F*F)~ 1 = 4l. Therefore (f>i = \<fi, for all % G J. In this case 
the reconstruction formulas above coincide and they look like 

16 J 
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for all / 6 Ji. In particular, moreover, if A = 1 and if all the vectors tpi are normalized, 
it follows that this frame forms an orthonormal basis of the Hilbert space, see |3|]. Of 
course, also the vice-versa holds true: if X is an orthonormal set in 7i, then X is a (1, 1)- 
frame of normalized vectors. This is an obvious consequence of the Parceval equality 
En | < <^n, / > I 2 = ll/ll 2 , which holds for all feH. 

How far can we go if the frame we are dealing with is not tight? In |3j and in the 
references therein, it is proposed a perturbative approach to this problem, which works 
well whenever A pa B. We will discuss in Section 4 a different approximation approach 
which can be used even when the frame bounds A and B are not very close. In the final 
part of this Section we summarize the standard procedure. We start defining a bounded 
self-adjoint operator R by 



R=t- F*F. (2.t 

A + B y 



Using inequality (|2.4|) we find that 



B-A B-A , , 

which implies that 

\\R\\ < (2.10) 

" " ~ B + A y ! 

Inverting definition (|278| ) we obtain F*F = ^±^(1— R), which also implies that = 
]4Xg(l — R)" 1 . Since certainly \\R\\ < 1 we can expand the operator (1 — -R) -1 into a 
norm-convergent series, Y^T=o R k - Therefore 

n 00 

^ = T-^E^Vr (2.H) 
A + ° k=0 

The N — th approximation is obtained considering only the first N — th contributions in 
the above infinite sum: 

2 N 

#=tt«5:*v ( 2 - 12 ) 

In P| it is shown that, for all / G TC, the following estimate holds: 



J6J 



< iiflir +i ii/n, (2.i3) 



which says that the above quantity converges to zero when iV — > 00. Of course, due to 
( p,10| ), this convergence is fast when A pa B, while it is rather slow if B — A is big, so that 
we cannot approximate (f>j keeping only few contributions in the expansion (|2.11|) . We 
need to consider iV big enough: this is the only way in which the error on the right hand 
side of ( |2.13| ) can be made as small as we like. A possibility to overcome this problem 
will be proposed in Section 4. 
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3 More Reconstruction Formulas 



We have seen in the previous Section that, given an (A, B)-frame X = {<p n }, it is possible 
to associate to it an unique dual (-^, -^)-frame X = {</5 n } and to use this two frames together 
to expand any vector of the Hilbert space. As it is well known, this is only one of infinite 
equivalent possibilities: in fact, since the sets X and X are not made up of independent 
vectors, in general, we cannot invoke any unicity of the coefficients of the expansion of 
a given vector. In other words: given a vector / G TC and a frame X, in general there 
exist infinitely many way to expand / in terms of (pj. Of course the possibilities given 
by Q2.6Q are particularly simple and elegant, so that they are considered of a particular 
interest. Moreover, these expansions solve the problem raised in the Introduction about 
the relation between the sequence {< <p n , f >} and the vector /, {(p n } being a frame. 

The problem we want to deal with in this Section is the following: does a given frame 
suggests other resolutions of the identity or, which is the same, other reconstruction 
formulas? 

The answer to this question is yes, and we will show how these different reconstruction 
formulas can be obtained and which techniques must be used. In particular we will show 
that the procedure discussed in the previous Section is only a particular case of this more 
general approach. 

Our first step consists in defining the operator 

Tx = F*F. (3.1) 

The norm of this operator is bounded from above and from below, A < < B, it is 

positive, T\ > 0, self adjoint, T\ = T\ , and its action on a given vector of Ti is given by 

Fif = EieJ < ( Pi,f > Pi- 

Let Ex be the family of spectral operators of T\. We can write, making use of the 

spectral theorem, 

< T x §^ >= [ B \d < Ex$,^ >, V$,#Gft. (3.2) 

J A 

Due to the fact that < A < B < oo, we can define arbitrary powers, positive and 
negative, of the operator T\. 

< J zr Q $, * >=< (F*F) Q $, * >= [ B \ a d < £ A $, * >, (3.3) 

J A 

for all and Va G R. 



6 



Following @ we see that T a , which of course still maps H into itself, is self adjoint 
and bounded for any a G R. Moreover we have: 

{FaY = J~a,p-) J~a.J~fi = ^a+/3- (3.4) 

Using the spectral theorem it is an easy exercise to compute the bounds for the operator 
T a . We find that 

K<\ <T 1 < 5 7 1 V 7 > (3.5) 
B^'t < ^ < A 7 1 V 7 < 0. (3.6) 

Given an arbitrary real number a let us define the following vectors: 

v\ a} = F a ¥i G J, (3.7) 

and let us call X^ the set of these vectors. We can prove that all these sets are frames 
in 7i. We have indeed the following 
Proposition 1. 

Each X^ is a frame. In particular X^ is an (A 2a+1 , B 2a+1 )-hame if a > — |, is a 
(1, l)-frame if a = — |, and is a (B 2a+1 , A 2a+1 )-frame if a < — \. 

Proof 

This result easily follows from inequalities ( |3.5|) and (|3.6|) and by the following equal- 
ities: 

E I < vt\ f > I 2 = E I < ^ f > I 2 = E I < > r = 

igj igJ iSJ 

= \\F(FJ)\\ 2 =< WJ,/> = < ^a+lfj >, 

which hold for all vectors / G 7i. □ 

Remark.- The above procedure teaches, in particular, how to get a tight frame with 
frame bound equal to 1 starting from a generic frame. Of course this does not imply that 
the vectors <yj- _1//2 ' ) form an ortho normal basis since normalization of these vectors is not 
ensured. 

We now define an operator 'associated' to the frame X^ a \ operator which plays 
the same role of F in the standard case. 

Let be an operator from Ti. into 1 2 (J) defined by the formula 

VfeH {F^f) j =<yf\f>. (3.8) 
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Since we have also < tpj, / >=< J- a (pj, f >=< tpj^Faf >= (F(J- a f))j, it follows that 

jritx) = F j7 a (39) 

Using now the equalities in ( |3.4| ) we see that, if we define the operator 

F {a) = F-x- a F*, (3.10) 

then 

for all a G R. These equalities can now be used to built up different reconstruction 
formulas, one for each value of a. First of all we observe that definition ( |3.8| ) implies that 

(^ (a) rc = £c^ (a) , (3.12) 
ieJ 



for all c G 1 2 (J). From ([3. 10 ) we have that, given / G H, 

(F {a) )*f = F(^ a f) = {< ipj^-i-af >} jeJ = 
= {< F.y^Pi, f >} jEJ = {< ^- l ~ a \f >} . £J • 

Therefore, given c G 1 2 (J), we have 

^c = Y,cM' 1 ~ a) - (3-13) 

ie J 

All these results can be collected to show that, for any / G 7i, we can write the following 
expansions: 

/ = £< V^J > ^' a) = £ < ^ l ' a \f > <pf\ (3-14) 

equalities which hold for all choices of the real a (Of course the second expansion is 
redundant since it can be obtained from the first one putting a —>■ — 1 — a). 

The proof is very easy: using formulas ( |3.8D , (|3.11|) and ( p,13|) we see that for any 



f eH, f = T^T^f = JF(«)(jr(«) = £. gJ < / > ^r l - a \ The other equality is 
proved in the same way, by considering the second equality in (|3.11|) . 



Let us now summarize our procedure: 

(a) the starting point is a given frame X = {y?-,} and the related frame operators F, 
F* and T x = F*F; 

(b) we construct the operator T a for a given real a as in (|3.3j ); 
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(c) we define a new frame whose vectors are (p[ a ^ = T a tpi 

(d) we also define its dual frame: <p\ 1 = T-i-^i] 



(e) we can now finally expand any / 6 H. as in (|3.14j) . 



Remarks. 

(a) It is interesting to notice that repeating twice the operation of taking the dual of 
a given frame we come back to the original frame. This result generalizes the analogous 
result for the standard procedure, tpj — (f>j, 0. More in details: if we take a frame 

and then we consider the new frame {f^} and its frame operator !F^ a \ it follows that 
the vectors constructed as in (£5|), {J^ a >J^)- x ipf \ coincide with ^~ 1_a) . In this sense 
we say that {fj -1 } is the dual set of {tpj}. Moreover, if we consider the "bi-dual" of 
{<fj}, then we obtain {f^} itself. 

(b) The method proposed here looks as a concrete version of a more abstract approach 
discussed in || where the attention was focused mainly on coherent states. 

(c) We observe that the approach we have followed here can be seen as a special case 
of a much more general one which is the following: given a frame {<fj} and its frame 
operators F and T\, considering the spectral decomposition ( |3.2| ), we can define new 
operators Fj = f(X) dE\, to be intended weakly, where / is any real bounded function 
with inverse bounded in [A, B\. In this way we could define a new frame \(p) } = {Ff^j} 
and work on analogous computations as the ones we have just discussed here, where we 
have always chosen /(A) = A Q , for some real a. 

Before discussing some examples of this procedure let us now consider some limiting 
case. We begin with the simplest one, a = 0. 

In this case we have (p[ a ^ = J^ofi = tpi and y9,-~ 1_Q ' ) = y?,- -1 '' = J--\<fi = <Pi, for all 
i G J. In this case the reconstruction formulas (|3.14j) collapse to the standard ones, (|2.6|) . 



Moreover the operators defined in this Section become: J-'^ = F, F^*' = (F*F) 1 F* 



F*, and so on. The equalities in (|3.11) coincide now with the analogous ones F*F 



F*F = 1, discussed in the previous Section in the standard setting. 

The situation is completely specular for a = — 1. In this case each (pf 1 ' coincides with 
y3j while y9 4 -~ 1_a ' 1 is nothing but y?j itself. Again the operators and T 1 ^ coincide with 
the ones widely discussed in the literature and in Section 2. 

Another limiting case is suggested by the Proposition 1 which assigns a particular 
importance to the choice a = — |. In this case we get a frame Z* -1 / 2 ) = {y?j -1 ^ 2 ' ) } which 
is tight with frame bound equal to 1. In fact, if a — — |, then — 1 — a = — |. For any 
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vector / G Tt the following equalities can be written 




(-1/2) 



and 



Ei<*r /2) ,/>i 2 = imi 2 - 



For what concerns the expressions of the operators we have, for instance, 1 ^ 
F ( F * F )(-i/2) and ^(-1/2) = ( F . F )(-i/2) F * = (jr(-i/2))*_ 



This last result is general: it is easy to show that for all a G R we have 




(3.15) 



Let us now consider some examples. Of course, the situation is much simpler for finite 
dimensional Hilbert spaces and it is just this situation that we start to consider. 

Example 1 

Let us take 7i = R 2 , and let us consider the three vectors ip± — (1, 0), ip2 — (0, 1) and 
ip 3 — -75(1, 1). It is easy to show that these vectors satisfy the following relation 



for all / G 7i, so that X is a (1, 2)-frame (the values of the bounds may be optimized, 
but we don't care about this point here). In the perspective of showing in which way 
our procedure works, it is convenient to start with frames which are not tight. This is 
because we just want to show how a (l,l)-tight frame can be obtained starting with an 
arbitrary other frame. It is obvious that J cannot be a Riesz basis, since its three vectors 
are necessarily linearly dependent. 



For any / = (f x , f y ) we get T x f = £f=i <<Pi,f>(Pi= l(3f x +f y , fx+Zfy)- Therefore, 
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||/|| 2 <EI<^/>| 2 <2||/|| 2 , 



i=i 



the operator T\ can be identified with the 2x2 matrix 
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where Ai = 1 and A2 = 2 are the eigenvalues of the E^s are the projection operators 
defined by E k g =< r} k ,g > k = 1,2, Vg G H, and 771 = ^-(1, -1) and r] 2 = ^(1, 1) 
are the eigenstates of JF X which correspond to Ai and A2. 

We are now ready to build up different pairs of dual frames, each for any given value 
of a. Let us start with a = — 1. 

We have T-\ — E\ + \E 2 and the set X^ 1 ^ = {j?-!^} is formed by the following 
vectors: y ? i~ 1 ' ) = |(3, — 1), <£>2 = \(~ 1;3) and <£>3 = 272^'^' ^ s we nave already 
discussed, is the dual frame of the original one X. In the standard language we 

would have written ip\ = (pi- It is straightforward to verify that, given / G 7~t, we 

o (—1) 3 (—1) 

can write / = Yh=i < V 9 ^ / > V 9 ) = X)i=i < V 9 ? >/ > V 9 *- Moreover we see that 
|ll/l| 2 < 5Zi=i I < fi\f > | 2 < ll/ll 2 ) estimate which does agree with the statement of 
Proposition 1. 

Of course the same results are obtained, mutatis mutandis, if we take a = 0. 
Let us now choose a = — \. 

The operator T_\ can be simply deduced by (|3.16|) : T_i — E\ + ^E 2 , so that 
we get the following vectors for the frame: (p[ 1 ' = |(1 + 4^, — 1 + 4^), (p 2 ^ = 
|(— 1 + 4=, 1 + 4g) and <£>3 = ^(1,1). It is a simple calculation to verify that these 
vectors form a (1, l)-frame. They are not, however, an ortho normal basis of 7i since none 
of these vectors is normalized. This, of course, is what we must have since the three 
vectors are necessarily linearly dependent in 7i. 

We conclude this example considering the case a = — |. 

The expression of the operator T_\ is again deduced by ( |3.16| ): T_\ = E 1 + ^3 E 2 , 



3 - ,< 3 4-, - 

and the vectors of the set T(-V3) are: = 1(1 + ^ -1 + <^~ 1/3) = |(-1 + 



7^73, 1 + -^3) and <y4 = ^7g(l, 1). The dual vectors are the ones corresponding to 
= _i + 1 = _|, that is ^" 2/3) = |(l+^3 , -I + W3), V [ 2 m = K-1 + W3 , 1+5W 
and (/?3 ^ = 5^(1,1). Again, it is straightforward to verify that the expansions in 
( ETni ) hold true, that is that, for any f e H, f = £f=i < ^ ( ~ 1/3) , / > <^ ( ~ 2/3) = £- =1 < 
V 9 ! f ~> V 9 ! 1//3 ^- Moreover, we can also check successfully the estimates given in 
Proposition 1: \\ff < £f =1 1 < J > ? < 2 1/3 ||/|| 2 and ^\\ff < £f =1 I < 

^ 2/3) ,/>| 2 < " 



Example 2 

We consider now another finite dimensional example a little bit more complicated. Let 
us take H = R 3 , and let us consider the four (linearly dependent) vectors <p\ = (1,0,0), 
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y?2 = (0, 1, 0), v?3 = (0, 0, 1) and y?4 = ^(1, 1, 1)- These vectors are a (1, 2)-frame: 

|2 <EI<^/>| 2 <2| ' 



i=l 



(Again, the frame bounds may be optimized). For any / = (f x , f y , f z ) G 7i we get now 

Tlf = im + fy + fz, f X + *fy + fz, f, + fy + ±fz), SO that 



/ 4 1 1 \ 



14 1 
V 1 1 4/ 

Using the spectral decomposition of T\ we can write 



Tx = E x + E 2 + 2E 3 



(3.17) 



where Ai = A 2 = 1 and A 3 = 2 are the eigenvalues of T]i = 4g(— 1,0, 1), r] 2 = 
77g(l, — 2, 1) and r] 3 = -^=(1, 1, 1) are the corresponding orthonormal eigenvectors and 
are the projectors on these vectors. We consider here only the case a = — |. 

In this case T_ i = Ei + E 2 + -^E 3 , and the (1, l)-frame is obtained as usual: ip\ 1 ^ 2 ^ = 



3 r _!(pi, for i = 1,2,3,4. We get <p 



(-1/2) 

i 



5(2+^,-1 + ^,-1 + 



V2> 



T2> 



<r2 



(-1/2) 



:- 1 + 75' 2 + 7I'- 1 + 7l)^3 = |(-l+^> -1 + ^,2+^) and** 



(-1/2) 



It is easy to verify that, for any / e Ti., f — Yn=i < >/ > V 9 

Et 1 |<^ ( - 1/2) ,/> 12 - 1 - 112 



,(-1/2) 



>/6 



1,1,1). 



and that 



Again, we observe that the various vectors are not normalized. Of course any attempt 
of normalizing these vectors would imply the loosing of their nature of (l,l)-tight- frame. 



Example 3 

We now consider an infinite dimensional example, previously discussed in 



Let 



us consider two real numbers po,qo G [ir, 2tt[. Let us then introduce A = and a 



C°°-function v{x) which is zero if x < and 1 if x > 1. In 
function 



it is discussed that the 



1/2 





sin 
1 

cos 



f^(+(- + x 

2 V A v po 



^ + A) 



/'u 



x < -—,x > — 

— po ' — po 

— < x < -— + A 

po — — po 

--+A<x<- 

Po — — PO 

— - A < x < — , 



I'll 



Po 



satisfies the following relations: 
( a ) Isupp g\ = 



12 



(b) Ekez\g(x-kq )\ 2 = ± 

Let now W(p, q) be the Weyl-Heisenberg operator defined by 

[W(p,q)f](x)=^f(x-q), 

for any square integrable function f(x). Given a function g(x) G 7i = C 2 (R) in reference 
||] it is shown that, by means of the above properties of the function g(x), the set {g mn } = 
{[W(mp , nq )g](x)}, for m and n both belonging to Z, is a tight frame: 2 mn6 z| < 
9mn, f > 1 2 = || /|| 2 for all f & TC. In this case it is particularly simple to write down 
the operator T\, which is simply T\ = -^1- Therefore, as it is obvious, the (1, l)-frame 
is formed by the vectors #<~ n 1/2) = T_^g mn = \p^9mn- 

Other non-tight frames can be obtained considering different values of a. 

A brief remark is in order: as it is clear from the above examples our procedure, which 
is quite general, can be easily implemented whenever we deal with finite dimensional 
Hilbert spaces or, when the original frame is tight. In both these situations, in fact, the 
operator T\ can be written in an easily handled form. On the contrary, for general frames 
in an infinite dimensional space, the situation, thought being clear from the theoretical 
point of view, is not as much easy to be implemented: the spectral projections must be 
replaced by the spectral measure and the finite sum which defines J-\ with an integral 
over a finite interval. T a is defined easily but, in general, it is not so clear which should 
be the explicit expressions of the vectors of the new frames, = T a <*Pi- 

There is a particular situation, which however is rather rare, in which the computations 
are easily done even in the infinite dimensional situation, namely when T\ is a finite-rank 
operator. In this case, in fact, the operator T\ can still be written as T\ = Y,k=i ^kEk, 
where is the k-th eigenvalue of T\ and E k the related projection operator. Of course, 
it would be interesting to obtain a condition on the original frame which ensures that 
the related operator T\ is of finite-rank. For instance it is clear that to any tight frame, 
and in particular to any orthonormal basis, it cannot correspond a finite-rank operator 
T\ since, being in this case T\ a multiple of the identity in 7i, its rank is necessarily the 
whole TC. 

In any case, this kind of problem is not new. It is just the same problem which we 
have to face with when we try to obtain T-\ = (F*F)~ 1 using the standard approach, 
0: except that in very few situations one has to proceed perturbatively, following, for 
instance, the procedure discussed in the Section 2. 

We will now show how this perturbation approach can be used for the computation of 
T a) at least when the frame bounds A and B do not differ too much. In the next Section 
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we will discuss a different perturbative approach which relaxes this condition. 

The problem we want to discuss is the way in which the dual operators T a and T-\- a 
can be computed perturbatively. We will restrict to some choices of a which have a 
particular interest. 

Let a G N. Using ( |2.8| ) we have T a = (1 ~ R) a i which is already given as 

a (finite) sum of positive powers of the operator R — 1 — -^-^F*F. More involved is 
the computation of T-\- a . Introducing the natural number M a = 1 + a, we see that 
J--i- a = (JF_ 1 ) a/ ° ! . Therefore, the problem of getting an expansion for T-x-o.-, is reduced 
to the analogous problem for which has been already discussed in Section 2. 

The situation is specular if we take a such that —a G N: in this case T-\- a can be 
found exactly, while the computation of T a is again reduced to the computation of T-\. 

Moreover, when a = 0, we go back to the standard procedure. 

For what concerns other values of a, we discuss here only the case of particular interest 
a = — j, which, as we have seen, gives rise to a (l,l)-tight-frame. 
Since \\R\\ < 1, we can expand the operator T_\ as 

Therefore we have, for any % G J, 

We define now an approximation of if { 2 of degree N: 

Given / G 7Y, how large is the error that we do when we substitute ^ 2 with 2 ' in 
the exact expansion / = < Pi 2 , f > 2 ? In other words, if we define the vector 

how big is the norm ||/ — /^||? This norm can be estimated telescoping it twice. Let us 
start introducing the self-adjoint operator 

T N = (l-R)-V*-p( ~1 /2 )(-R) k . 
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Therefore, 



1 T^ifi = tp i 2 — <fi 2 , for all i E J. It is easy to see that 



A+B 



ll/-/ (JV) ll<IIE<(^ H) -^ H),iV U>^ H) ll + 



ieJ 



+ 



+llg<(^ ) ' N -^ ) ),/X^' ) -^ ) ' N ) 

+ ilE<^» l ,/x^» , -^ ) '")il, 

which shows that the operator T/v plays a relevant role in the estimate of this norm. 

Let us consider in some details, for instance, the computation of the first contribution 

above. We first notice that < ((p\ ~ 2) - tpf 1 * ),f >= ^JjTb < T NPh f >= \Ja+b < 

(-1) 

(pi, T/v/ >• Moreover, using the notation of Section 2, we also deduce that < g, 2 >=< 
g, T\ii(pi >=< J-'i^g, <fi >■ Therefore, this contribution can be estimated in the following 
fashion: 



/ (-3) (~h)> N \ r (-5)11 

E<(^ )J>Pi II 



ieJ 



A + B m<i 



sup 

IMI<i 



A + B' 



A + B 



E < Pi,T N f >< Tx/ 2 g,Cpi > 



Fl/2T N f\\ < 



2B 



A + B' 



IT 



N 



Here we have used, among the other things, the fact that: ||T\/ 2 1| < yB, which directly 
follows from ( |3.5|) . The other contributions above can be estimated in analogous way. 
The final result is the following: 



\\f-f W \\< 24 



25 2B llm . 

+ -r-^\\T N \ 



IT 



A' 



A + B A + B' 

We see here the necessity of estimating ||T/v||. The details of this computation are given 



in the Appendix, and the result is in formula (|A.2|) . Here we give only the result: 

< B _ A ^N+l " 



Tv < 



2A 



I A + B 
2A ' 



The conclusion is therefore the following inequality: 



Id / d A\ N+l 

\\f-f iN) \\<\R r 



A V 2A 



IB {B-A- N+l 
2 + 4,/ — 



A \ 2A 



which converges to zero in the same fashion of the analogous estimates given in ||, at 
least if B < 3 A. Obviously the speed of this convergence increases as much as A becomes 
closer to B. 
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We want to end this long Section with some reversibility remarks. The first question 
is the following: given a frame, it gives rise to a single operator T\ = F*F. Is this 
correspondence a bijection? That is, given a self-adjoint operator T\ which satisfies the 
properties discussed at the beginning of this Section, is it possible to associate to this 
operator a single frame? 

It is easy to find counterexamples which show that in general this is not so. The 
first trivial case is the following: if the set {y?i} is a frame then also {e toti (fi} is a frame 
for any choice of the real constants aj, and both the frames correspond to the same 
T\ operator. Therefore it is not possible, given !F\, to fix a single 'generating frame'. 
However, it could happen, in principle, that T\ fixes only a frame of rays in the Hilbert 
space, and not a frame of vectors. This cannot happen, in general, for finite dimensional 
Hilbert spaces as we can show with a very easy example: let us suppose that dim(7i) = 2 
and that TC is a real Hilbert space. Then T\ is a 2 x 2 real and symmetric matrix with 
three independent entries. Let us now consider a frame of TC of four vectors. Condition 
T\j — J2i=i < <Pi, f > <Pi f° r an y / ^ gives three independent equations which become 
seven if we also consider the normalization conditions of the vectors. They are not enough 
to fix all the vectors of the frame! 

The question is still open for infinite dimensional Hilbert spaces. 

The second remark is very much related to the previous one: given an (A, 5)-frame 
{(fi} and its related operator T\, we find an unique (1, l)-frame {^p'f 1 ^ 2 ^}- The question 
is again a problem of reversibility: given now an (1, l)-frame {<pi}, is it possible to find an 
unique (A, 5)-frame {^} such that & = ^ { ~ 1/2) ? A gain, the answer is in general negative 
and it can be shown simply by giving examples of (1, l)-frames each of which can be 
obtained by different non-tight frames. Before discussing these examples we remind that 
normalization of the vectors of the frame is very important: it is enough to think to the 
examples of (1, l)-frames discussed in this Section, which are certainly not ortho normal 
bases due to the lack of normalization of their vectors. It is well known, in fact, that such 
a frame turns out to be an orthonormal basis if its vectors are all normalized, || . 

We start observing the following property: 

- If {y?i} is a frame then also = oupi} is a frame and we have (p\ 1//2 ' ) = ^ for 
alii e J and Va > 0. 

The proof of this statement is a straightforward computation. With obvious notation 
we call _F$ and F v the frame operators related to the two different frames. By the definition 
we have, for any / e TC: 
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so that = a 2 F*F v . Therefore we have 

This is a particular case of a more general situation which is easily considered if dim 
7i < oo. Let {tpi} be a given frame, with % — 1, ...,N, and let i 7 ^ be its frame operator. 
Let us suppose that it exists a self-adjoint operator A e B(7i) which commutes with 
F*F V and for which there exist two positive real numbers < A_ < A + < oo such that 
A-TL < A < A + H. In general any operator of the form al satisfies these conditions for 
any strictly positive a, independently of the form of J 7 ^, and this is the case we have 
considered above, as we will see in a moment. In general such an operator, if it exists, 
must be related to the frame operator. Let us now define $; = A 1 / 2 (p i for all i = 1, N. 
The set of these new vectors is a frame whose related frame operator is F$ = F^A 1 ^ 2 . We 
have: 

$r /2) = (f^- 1 ^ = (F^y-^A-^A 1 ^ = (w ( - 1/2 Vi = v4 _1/2 \ 

as it was to be shown. The extension to the case dim Ti = oo contains some subtle points 
and it will not be considered here. 



4 A Perturbation Approach 

We have discussed in the previous Sections the difficulties of finding the operator T-\ = 
(F*F)~ 1 for a given frame, in general. The way in which this problem is overcame in the 
literature is to consider a perturbative approach whose speed of convergence, however, 
depends strongly on the difference of the frame bounds, B — A: the smaller this quantity 
the more rapid the convergence of the procedure. In this Section we will propose a different 
perturbation approach which relax this condition, so that it can be also used when B — A 
is not so small. 

The idea is very simple and it is strongly related to the fact that the operator T\ = F*F 
is strictly positive. Let {</?j} be a frame and F its associated frame operator. We know 
that 

At < F*F = Tx < Bl. (4.1) 

In all this Section we will not assume that is much smaller than I. On the contrary, 
the procedure we are going to discuss is particularly useful in the case in which this ratio 
approaches the unity, that is in the case in which B is much larger than A. How can we 
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construct a perturbative approach to get the operator T-\ = (F*F) x ? We will propose 
a procedure which depends on the values of the frame bounds. 



Let us first assume that 

1 < A < B < oo. (4.2) 

As it was discussed in the previous Section, due to the nature of the operator T\ it is 
possible to find a spectral decomposition for this operator and we can build up functions 
of this operator. For instance, recalling equation ( |3.1| ) 

< >= [ B \d < £ A $,^ >, V$,$6H, (4.3) 

J A 

we see that the operator Log B (jF 1 ) is defined by 

< Log B (^ r i)$, \& >= f B Log B (A) d < E\&, \1/ >, V$,^Gft. (4.4) 

J A 

The integral above is surely well behaved since A > A > 1. Using functional calculus it is 
easy to verify that inequality ( f4. 1|) gives rise to 

at < Log B (^i) < 1. (4.5) 

This result is, of course, a consequence of the fact that AH and BTL both commute with 
T\. In ( fL5| ) we have defined a = Log B (A). Of course < a < 1. The relevant difference 
with respect to the original double inequality ( |4.1|) is that the two new bounds a and 1, 
are necessarily near to each other. We define, in analogy with the standard approach, an 
operator R via 

fl = l_-J-Log fl (7i). (4.6) 
1 + a 

Using the bounds (|4.5|) we find that the operator R satisfies the following double inequal- 
ity: 

1 a i<i*<i^lL. (4.7) 



1 + a 1 + a 

Equation (|4.6|) can now be inverted and it gives, after some minor calculation, 

T_ x = (F*F)~ 1 = -^B^ R = —^e laBl ¥ R . (4.8) 

Vab Vab 

Let us now see how the perturbative approach works. We call, as usual, (f>j = F^iipj, 
which is the exact dual frame, and we define its N-th approximation by 
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We are interested in computing which error do we do if, instead of < V 9 ^ / > &i = /? 
we compute = Y^i^j < fit f > 0f- m other words, we want to estimate ||/ — f^\\. 
The procedure is similar to the one we have discussed in the previous Section. Let us 
define the self-adjoint operator Zn, 



Z N ^e^ R -j:^(lnB 1 -±^R 



Using the definitions and the following equality 

1 



~ N 

< <Pj -<Pj ,g > 



sfAB 

where, as usual, <f>j = T-ifj, we have 

\\f-fW\\= sup 

llsll<i 



VAB 



sup |< f,T\Z N g >| < J— 



Il3ll<l 



.4 



The problem is, therefore, to evaluate ||^jv||. In the Appendix we will deduce that 

1 



Z N \\ < 



A (N + l) 



I- a n N+l 
- — -lnfl 



so that 



N+l 



11/ - /'"'II < 1 11/11 ('f^, . 



(4.10) 



which shows that ||/ — f^\\ is exponentially convergent to zero when N — > oo, indepen- 
dently of the values of A and B, which can also be very different from each other. 



Let us now assume that 

< A < B < 1. 

In this case the inequality ( |4.1| ) can be rewritten as 

01 < Log^) < 1, 



(4.11) 



(4.12) 



where f3 = Log A B is strictly positive and (3 < 1. We define the operator R in a way which 
is similar to what we have done in (|4.6|) : 

2 



R = l 



1-/3 



Log A (^"i). 



(4.13) 
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The operator R is bounded, ||i?|| < ^—1, and the operator T-\ differs from the one given 



1+/3 

in (14.81) and turns out to be 



= ^=A^ R = _^e lnAi ^ fi . (4.14) 

Again, if we define in the usual way, we get for ||/ — f^\\ an analogous estimate 
as the one we have obtained before, showing its exponential convergence to zero for iV 
diverging. 

The last condition on the frame bounds A and B is 

0<j4<1, B>1, (4.15) 



which can be reconduced to the first situation simply by multiplying ( |4.15| ) for, say, -§■ 
In this way this inequality can be rewritten as 

21 < jFi < bl, (4.16) 



where b = At this point we call 5 = Log b 2, so that Q4.16j ) becomes 



51 < Log b (^i) < 1, (4.17) 

and we can proceed exactly as we have done before. We obtain, after some minor com- 
putation, 

= -^b^ R = -]=e l « bl -¥ R . (4.18) 

The estimate of \\f - f {N) \\ goes through as usual. Again, this quantity converges to zero 
exponentially with N. 

The following remark is in order: while in the standard approach the zero approxima- 
tion of (fj, <p®, is nothing but with this approach the zero-th approximation turns 
out to be -7= = (fj . The arithmetic mean is replaced with the geometric one. 

Of course the standard perturbative approach we have used in the previous Section in 
order to discuss how to obtain could be replaced with the different approach we have 
just developed. However, we will not consider this analysis here. 



5 Frames Versus Orthonormal Bases: Some Remarks 

This Section, which concludes the paper, is devoted to discuss some differences between 
a set of orthonormal vectors and a frame. 
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As we know, [0, a set {<fij}, j G J, is an orthonormal basis of a given Hilbert space 
7i, if and only if one of the following equivalent conditions is satisfied: 

(i) for all / G U then ||/|| 2 = E ie j I < f,(pj > | 2 ; 

(ii) for all / G H then / = £ ieJ <<pj,f> <Pj m , 

(iii) for all f,g eH then < f,g>= E je j < f, <Pj >< <Pj, 9 >; 

(iv) if / ETC is such that < /, <£>j >= for all j G J, then / = 0. 

Similar conditions can also be restated for an (A, S)-frame {0j} with frame operator 
F . Introducing the dual frame of vectors 4>j = (F*F)~ 1 <ftj, the above conditions become: 

(i') for ail feU then A||/|| 2 < £ ie j | < /, <j>j > | 2 < S||/|| 2 ; 

(ii') for all / G n then / = £ jeJ <(f) j ,f> 4>f, 

(iii') for all f,geH then < f,g>= EjeJ < /> 0j >< 0j; 9 >5 

(iv') if / G W is such that < /, 0j >= for all j G J, then / = 0. 

These conditions can also be rewritten in terms of the frame operator F. In this way 
conditions ii') and hi') collapse and the conditions are: 

(a) Al < F*F < Bl; 

(b) the operator (F*F) can be inverted and its inverse (F*F)~ l is bounded in 7i; 

(c) ker F = {0}. 

The equivalence between (a) and (i') is well established in the literature. For what 
concerns the equivalence of condition (b) with the conditions (ii') and (iii'), this is easy 
to be deduced: since the operator (F*F) _1 belongs to B{7i), we can define a new op- 
erator F = F(F*F)~ 1 which maps H into 1 2 (J) such that F*F = F*F = 1. Defining 
(f)j = (F*F)~ 1 (pj one direction of the equivalence is proved. The converse implication is 
immediate, since the existence of the dual frame already implies the existence of (F*F)~ 1 
and its boundedness. Finally, if ker F = {0}, then the operator F is injective, so that 
Ff = only if / = 0. Vice-versa, if Ff = implies that / = 0, then F is injective as an 
operator from n into 1 2 (J), so that its kernel contains only the zero vector. 

It is widely discussed in the literature that condition (a) implies condition (b), or, in 
other terms, that (i') implies (ii') and (iii'), essentially as a consequence of the inequality 
F*F > AIL. It is also very easy to see that condition (i') implies condition (iv') or that (a) 
implies (c): let us suppose that / G n is orthogonal to all the frame vectors. If condition 
(i') holds it is obvious that ||/|| = 0, so that / = 0. 

The main point is now to show that, contrarily to what happens for orthonormal bases, 
the converse implications are not all true. In particular, we will show that condition (b) 
implies condition (a), while condition (c) is not equivalent to the other two conditions. 
In order to show the implication (b) =^> (a) we observe first that the operator (F*F)~ 1 is 
positive as its inverse. Therefore, using the fact that if the bounded operator T is positive 
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then T\\T\\ > T 2 , ||, we deduce that the operator F*F satisfies the following inequalities: 
H( F ,^-i|l 1 < F*F < Therefore F is a frame operator. 

The fact that the last condition does not imply, say, the first one, follows from the 
fact that the condition of F being injective cannot of course imply also its surjectivity. 
Therefore the set of vectors related to this F in general does not generate the whole 
Hilbert space, so that it is not a frame, in general. On the other hand, if we also require 
to F of being surjective, then it would be invertible, condition which we know is not 
required at all to any frame operator. This is what happens for Riesz bases which, as 
we have discussed in the Introduction, generate the whole Hilbert space and contain only 
linearly independent vectors, 0]. 

We conclude, therefore, noticing that, even if the conditions which define an orthonor- 
mal basis can be easily extended to general frames, these new conditions are not in general 
all equivalent. This difference is essentially due to the linear dependence of the vectors of 
the frame. 
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Appendix : Norm Estimates 

This Appendix is devoted to show the details of the estimate of the norms of a couple of 
operators which are used in the paper. 
Let us start with the first operator, 



which appears in Section 3 in the computation of ||/ — /^||- The estimate proceeds 
using the functional calculus in the following way: let us introduce the spectral family E\ 
of the bounded and self-adjoint operator R. Of course this spectral family coincides with 
the one of the operator T\ = F*F, see equality ( |2.8|) . As in equation ( |3.2|) , we can write 





(A.l) 
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where 7 = ^tj- Being (1 — R) a positive operator bounded from below from a strictly 
positive constant, we can consider the self adjoint operator (1 — i?)~ 1//2 , 0, which can be 
written, in terms of the spectral family E\, as 

< (1 - i?)~ 1/2 $, >= r (1 - \)- 1/2 d < £ A $, >, V$, *GH. 

Obviously we have also 

< E ( -1 / 2 ) i-R) k ^>= [_ E ( _1 / 2 ) (-A) fc rf < s A $,« > . 



Therefore, for any / G H, we get 

2 

d<EJ,f>. 



\T N f\\ 2 ' 



;i-a)- i/2 -e( ! /2 K-a) 



A 1 

E 

fc=0 

The problem is now reduced to estimate the following function 



k 



p(A) = (1 - A)~ 1/2 - £ f 1/2 ] (-A)*. 



The idea is to observe that p(A) is nothing but the N-th (Lagrange) remainder of the 
Maclaurin approximation of the function ^(A) = (1 — A) -1 / 2 , so that we can estimate 
p(A) in the usual way: p(A) = ^v+ffi X N+1 . Here £ is a suitable real belonging to [0, A] if 
A > or to [A, 0] if A < 0. 

2JV+1 

With some easy estimates we see that ty( N \£) < N\ (jzr) 2 , which implies that 

2JV+3 

|p(A)| < (^-) 2 ^ N+1 \ Therefore we get 



27V+3 / 1 \ 27V+3 



( 1 ,7 / 1 V 



7 



and, finally, 



The second quantity we want to estimate in this Appendix is the norm of the operator 
Zjv introduced in Section 4, 
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The approach is exactly the same as for the previous norm: we use functional calculus, 
Maclaurin expansion and the Lagrange formula for the remainder. We do not give all the 
details here, but only the result which is the following: 



>B±J2||fl|| 

(N+l)!" V 2 ~* J -(N + 1)\\A 

(A.3) 

where a = Log B A and, due to Q, \\R\\ < ±=%. 
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